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INTRODUCTION 

Theoretical  models  ot  the  D-region  have  increasingly  relied  upon  com- 
plicated computer  code  simulations  of  the  chemical  kinetics  and  photo- 
chemistry of  the  region.  Recent  progress  in  developing  l)-region  models 
and  some  of  the  problems  inherent  in  verification  of  the  existing 
models  have  been  outlined  by  Sechrist  [1]  and  Aikin  [2].  As  both 
authors  point  out,  progression  of  the  complexity  of  D-region  models 
from  the  traditional  lumped  parameter  models  iwhich  require  only  single 
positive  and  negative  ionic  species  plus  electrons)  to  the  multispecies 
codes  places  an  almost  impossible  task  upon  the  experimenter  who  seeks 
to  verify  such  codes.  Nevertheless,  the  failure  of  lumped  parameter 
models  to  adequately  describe  the  behavior  of  the  D-region  under  even 
the  most  controlled  conditions  imposes  the  requirement  for  development 
and  verification  of  more  physically  realistic  kinetic  chemistry  models. 

This  report  presents  the  results  of  a studv  which  was  undertaken  to 
partially  verify  the  D-region  chemistry  code  known  as  DAIRCHEM  (acronym 
for  D-region  air  chemistry)  which  models  the  effects  of  time-varying, 
naturally  occurring  ionization  sources  on  the  chemistry  of  a parcel  of 
air  in  the  middle  atmosphere.  The  emphasis  is  on  the  ion  chemistry 
with  sufficient  neutral  chemistry  to  make  a self-consistent  model. 

Aikin  (2]  has  pointed  out  that  verification  of  D-region  models  requires 
data  on  not  only  the  spatial  and  temporal  distributions  of  the  ambient 
ion  and  neutral  species,  but  also  on  the  Ionizing  source  functions  as 
well.  The  data  base  used  in  this  study  consists  of  an  assembly  of 
ionizing  source  data  and  models  and  I'-region  data  obtained  from  a wide 
variety  of-  experimental  measurements  made  during  the  12  November  1966 
eclipse.  A preliminary  study  which  gathered  these  data,  placed  them 
in  a common  coordinate  svstem,  and  defined  the  ionizing  source  func- 
tions, was  conducted  by  Sears  [3], 

Comparisons  with  experimental  data  provide  a means  of  validating  the 
basic  modeling  computer  codes  which  are  used  as  input  for  Army  communi- 
cation svstems  and  In  the  Army  nuclear  weapons  effect  community.  In 
addition,  such  comparisons  are  a means  of  assessing  what  further  ex- 
perimental measuiements  must  he  made  to  provide  data  to  improve  already 
exist ing  computer  codes. 

The  organization  of  this  report  is  as  follows:  "Ionospheric  Measure- 
ments" describes  in  detail  the  1966  solar  eclipse  experimental  measure- 
ments of  positive  and  negative  ion  densities  versus  altitude  and  time. 
"D-Reglon  Ionizing  Sources"  defines  the  ionizing  source  functions 
appropriate  to  the  circumstances  of  the  eclipse,  and  "Description  of 
the  DAIRCHEM  Computer  Code"  briefly  describes  the  DAIRCHEM  computer 
code.  Under  "Comparisons  of  DAIROIEM  with  the  1966  Eclipse  Data 
Base...,"  the  Ionizing  source  functions  are  incorporated  into  the 
DAIRCHEM  code  and  the  predicted  ion  and  electron  densities  are  com- 
pared with  the  experimental  data. 
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IONOSPHERIC  M E AS  U REM  ENTS 
D.it.i  Sources 

Many  ionospheric  measurements  were  made  in  Brazil  and  in  other  South 
American  countries  on  tlu  path  of  totality  of  the  1966  eclipse.  The 
larger  portion  of  these  measurements,  however,  were  made  individually 
and  not  in  coordination  with  other  ground,  aircraft,  or  rocket 
launched  experiments,  and  hence  are  of  only  nominal  value  to  the  pres- 
ent analysis.  The  Defense  Nuclear  Agency  (DNA)  (formerly  Defense 
Atomic  Support  Agency)  supported  experiments,  plus  those  of  National 
Aeronautics  and  Space  Administration  (NASA)  and  local  Biazilian  experi- 
menters which  were  carried  out  in  the  vicinity  of  the  Rio  Orande  rocket 
launch  complex,  were  closely  coordinated  with  one  another  and  wit  1 pro- 
vide the  major  fraction  of  the  data  which  will  he  described  and 
summarized  in  this  section.  The  experiments  which  are  considered  in 
detail  here  are  summarized  in  tables  1 and  J. 

Three  classes  of  data  exist:  ground-based  and  a ireraf t -based  measure- 
ments of  electron  density  versus  altitude  which  are  derived  from 
analyses  of  propagation  effects;  in-situ  measurements  of  electron 
density,  positive  ion  density  and  species  composition,  and  electron 
temperature  which  were  carried  on  rocket-launched  probes;  and  roeket- 
to-ground  propagation  experiments  which  yielded  altitude  profiles  of 
electron  density.  In  this  portion,  the  authors  take  advantage  ot  the 
redundancy  available  in  the  electron  density  data  to  attempt  to  provide 
a best-fit  to  the  time  and  altitude  variations  of  electron  densitv 
during  the  eclipse.  Because  different  measurement  techniques  have 
different  analytical  and  empirical  sources  ot  error,  a combination  ot 
two  or  more  sources  of  data  should  in  principle  yield  results  which 
are  more  truly  representative  of  the  D-region  during  the  eclipse. 

Eclipse  Coordinate  System 

The  principal  difficulty  in  combining  data  from  ground,  aircraft,  and 
rocket-based  measurement  platforms  is  the  accurate  specification  ot 
their  location  with  respect  to  a common  eclipse  oriented  coordinate 
system.  To  provide  the  common  coordinate  system  required,  a system 
utilizing  the  fewest  number  of  independent  variables  possible  was 
chosen.  The  significant  variables  are  altitude,  eclipse  t ime,  and  the 
duration  of  totality  at  the  chosen  altitude.  Eclipse  t ime  is  specified 
as  seconds  before  second  contact  or  seconds  after  third  contact  for 
measurements  outside  totality,  or  as  seconds  after  second  contact  and 
duration  of  totality  for  measurements  within  totality. 

The  locations  of  the  trajectories  of  several  DNA  and  NASA  sponsored 
rockets  within  the  altitude-time  coordinate  system  are  plotted  in 
fig.  1.  Trajectory  information  for  the  NASA  rockets  in  this  coordi- 
nate system  was  taken  from  Mechtly  et  al.  [4].  This  coordinate  system 
allows  one  to  avoid  the  problem  of  location  of  the  data  in  horizontal 
coordinates  at  a given  altitude  within  the  eclipse  totality  because  the 
pertinent  coordinates  are  time  after  second  contact  and  the  duration  of 
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TABLE  1.  ROCKET  PROBE  MEASURi^ENTS  OF  ELECTRON  AND  ION  DENSITY 
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totality.  Sears  [3]  gives  a more  detailed  account  of  placing  the  tra- 
jectories within  the  eclipse  coordinate  system. 

Rocket  Probe  Fleet ron  Density  Measurements 

Measurements  of  electron  density  were  made  by  a wide  variety  of  rocket- 
borne,  in-situ  probes  throughout  the  eclipse.  In  addition,  two  certi- 
fication tilings  were  conducted  on  5 November  I9bb,  a week  before  the 
eclipse,  which  obtained  background  data  at  the  same  time  of  day  as  the 
eclipse.  In  all.  Id  successful  firings  were  made  during  the  eclipse 
period.  The  electron  density  instrumentation  capabilities  of  these 
rockets  and  their  launch  times  are  summarized  in  table  1.  In  this 
section,  tiie  measurements  ot  electron  density  are  summarized  and  placed 
into  the  eclipse  coordinate  system  previously  described  so  that  tem- 
poral comparisons  mav  be  made  of  electron  density  at  a given  altitude. 

Ttie  different  experimental  techniques  used  to  measure  electron  density 
from  a rocket  vehicle  during  tiie  19t>6  eclipse  include:  Langmuir  probes, 
subsonic  blunt  probes,  and  rocket-to-ground  propagation  measurements 
to  calibrate  the  Langmuir  probe  measurements  [*'»].  Other  experimenters 
used  absolute  laboratory  calibrations  for  probe  measurements  [Si.  In 
each  case,  conditions  mav  exist  wherein  the  calibrations  adopted  for 
the  Langmuir  probes  during  the  eclipse  mav  be  erroneous.  In  addition, 
propagation  experiments  which  are  analyzed  on  the  basis  of  a strictly 
horizontally  stratified  ionosphere,  ignoring  horizontal  or  temporal 
variations  in  electron  density,  mav  also  be  in  error,  Ir.  the  case  of 
noneclipse  conditions,  most  of  the  measurement  techniques  produced 
closely  similar  measurements  of  electron  density  versus  altitude  with- 
in the  probable  errors  in  their  respective  calibrations.  Figure  2 
illustrates  the  background  electron  density  versus  altitude  curves 
obtained  bv  the  variety  of  rocketborne  probes  summarized  ir.  table  1. 

As  shown  in  tiie  figure,  most  measurements  above  about  70  km  agree 
within  20  percent  or  so  for  the  background  conditions.  Below  70  km, 
blunt  probe  data  appear  to  diverge  from  the  other  results.  This  dis- 
crepancy has  not  been  investigated  in  detail.  To  emphasize  the  posi- 
tive aspects  of  this  comparison  of  various  experimental  techniques, 
the  NASA  Langmuir  probe  data  (calibrated  by  propagation  data)  , the 
Ballistic  Research  Laboratories  (BRL)  low  frequency  (,LF)  propagation 
experiment,  and  the  BRL  Langmuir  probe  measurements  of  electron  den- 
sity all  agree  above  about  t>5  km  within  their  experimental  errors. 

Below  P5  km,  where  the  electron  density  is  smaller  than  100  cm--,  the 
individual  experimental  errors  are  as  large  as  a factor  of  3 or  so, 
thus  estimates  of  the  absolute  electron  density  curves  below  this 
level  may  possess  relatively  large  experimental  error. 

The  most  significant  data  which  were  derived  by  this  study,  after 
putting  all  of  the  rocket  probe  data  on  a common  eclipse  coordinate 
system,  are  illustrated  in  fig.  3.  The  electron  density  versus  time 
measurements  at  selected  altitudes  between  70  and  90  km  are  depicted, 
utilizing  all  of  the  available  rocket  probe  and  rocket-to-ground 
propagation  data.  Three  significant  features  are  to  be  noted  in  this 
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BACKGROUND  ELECTRON  DENSITY  PROFILE 


Figure  2.  Background  measurements  of  electron  density  versus  altitude  from 
in-si tu  rocket  probes „ LP* Langmuir  probe,  BP=Blunt  probe,  and 
P=*rocket- to-ground  propagation  experiments 
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Figure  3.  Electron  density  versus  time  about  totality. 


ilputo:  the  relative  rate  ot  clianpo  aiul  mapnltudc  ot  ilumia'  in  eleo- 

t ron  ilous  1 1 v I'otwoon  t list  anil  noat  sooonil  contact  ts  small,  tin?  rato 
ot  ilooroaso  In  electron  aonsllv  immediately  following  sooond  oontaot  is 
votv  large,  ami  t ho  absolute  electron  ilonsltv  ts  reduced  within  JO 
sooomls  1 1'  mi  tint  oh  hv  touphlv  a laotor  ol  A to  10  in  contrast  to  tho 
gradual  ilooroaso  oxhlhitoil  Online  tlto  initial  stages  ot  tho  eclipse. 

To  hi  ln>',  out  this  strong  totality  related  ioatuto,  which  was  suggested 
hv  tlio  NASA  locket  data  published  hv  Moohtlv  ot  al.  [ •<  ] , all  ot  tho 
rocket  pi  oho  data  wore  used,  and  tho  o loot ron  density  values  were  taken 
as  tho  ptinclpal  expel  intent  c i s insults  indicated. 

The  impel  taut  ionesphoite  hehavioi  which  has  hceu  explicitly  demonstrated 
in  tills  section  ts  the  strong  loduotion  in  elootron  density  inimed  i a t e 1 v 
around  tho  togion  ot  totality.  1 yen  it  the  int ei oompat Isons  ot  the 
various  rocket  probe  results  illustrated  hv  tip,.  ' aio  do  emphasized  hv 
asst  put  up  a laipo  prohahle  eiioi  to  each  pi  oho  data  point,  t lie  chaise 
teristle  rapid  decrease  in  electron  density  at  second  contact  is  pto- 
soived.  Ionospheric  data  I row  othoi  eclipses  have  also  suppested  this 
ottect  | (' , ' ] in  the  70  to  'hi  km  altitude  i opinio.  Tho  implications  of 
these  data  with  respect  to  atmospheric  do i on ( rat i on  model i up  ai o dis- 
cussed mulct  "I'ompa i t sons  ot  1'Al  Ki'IU'.M  with  the  l‘h>t>  |,  lipst  Oats..." 

To s i t i v e 1 on  Pens i t \ Pa t a 

Positive  (on  density  and  species  vote  measured  b\  in-situ  locket  probes 
over  tin'  altitude  ranpo  of  M'  to  100  km.  Three  types  ot  probes  weie 
used:  a parachute  deployed  blunt  probe  to  measure  total  posit ive  ion 

density  below  about  SO  km  [8],  a ret  aril  (tip  potential  analyzer  | 0 1 , and 
a positive  ion  mass  spec  t rottiet  ei  to  measure  total  positive  ion  density 
ami  ionic  species  mass  distribution  above  about  '0  km  [ 10,11).  A 
small  common  altitude  taupe  tabeut  0 to  80  kirO  was  .nailable  to  both 
the  blunt  probe  and  the  ion  mass  spec t l emit ei . 

As  desotibod  hv  liakot  | 8 | , tho  positive  i oiuluc  t i v i t v et  tho  ionospheio 
is  determined  t rom  the  euiieut  to  the  nop.it  tvolv  charped  blunt  probe. 

The  positive  ion  niimbei  density  is  iclated  to  the  measured  conductivity 
o bv: 


who  i o 

N is  tho  positive  Ion  density 
e is  tho  electron  eharpe 

i,  is  the  reduced  mob  i 1 1 1 v of  the  positive  ton. 

The  mobility  value  chosen  tot  l nt erpt o t at l on  of  the  data  was  n - 1.8 
cm'V'hi  , which  is  typical  ot  values  piovlouslv  moasui ed  t oi  the  0, 
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ii'ii.  llowovoi  , us  Nan  (si's  mass  spoo 1 romot  ot  mousut  omont  s showml , tho 
pi  I m l pa  I posit ivo  Ion  is  II  yO*  i II  ,0)  , im  hi  phot  hviliatos,  Is  low  SO  km  .u 

so,  tlio  nii'l'  I I 1 1 v v'l  wlttoli  is  takou  as  ,’.S.  Whou  this  inoroasoil  moh  i I 
Itv  valuo  is  iisoil,  tin  posit  iw  ion  ilonsitios  ropoitoil  t 01  t ho  Munt 
I'toho  ospoi  imont  aro  somowhat  toiluooil.  A summary  ot  t ho  ail  fust  oil  hlunt 
pioho  ilata  ami  t ho  oomhinoil  totaiilinp  potontlal  .inalvzor  (.Kl'M  am1  ion 
mass  spootromotoi  total  ion  Sonsitv  ilata  is  oontaino.l  in  tip.  i.  Ms  vo 
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Figaro  t>  shows  a comparison  ot  the  partial  reflection  electron  density 
data  with  in-sltu  rocket  probe  measurements  at  several  times  during 
the  eclipse.  These  curves  show  that  the  partial  reflection  tesults 
agree  closely  with  the  rocket  borne  Langmuir  probe  data  obtained  bv 
HKL,  Air  Force  Geophysics  laboiatorv  (formerly  Air  Force  Cambridge 
Research  Laboratory),  and  NASA  instruments. 

The  good  quantitative  agreement  of  partial  tef lection  sounder  and 
rocket  probes  utilizing  a variety  of  experimental  techniques  must  be 
considered  as  verification  of  the  accuracy  ot  each  technique. 

examination  ot  the  ground-based  electron  density  data  indicates  that 
such  experiments  have  two  valuable  roles:  cross-calibration  and  en- 
hancement ot  rocket-based  N profiles  versus  altitude  and  time  bv 

direct  sounder  measurements;  and  provision  of  a basis  for  comparison 
of  eclipse  O-region  qualitative  phenomenology  with  prior  efforts  using 
indirect  techniques  lor  Inferring  O-region  electron  densities,  e.g., 
riometers  and  VLF/l.F  sounders. 

0-KKG10N  ION  ICING  SOURCES 

Maintenance  ot  t he  ionization  in  tho  duvllt  O-roglon  between  about  Ml 
and  *>()  km  requires  several  sources  of  ionization.  Tho  intensity  of 
each  source  has  its  particular  dependence  upon  temporal,  geographical, 
solar  cvcle,  altitude,  and  other  variables  each  ot  which  must  be  quan- 
titatively specified  for  the  circumstances  of  the  l‘lbn  eclipse.  Addi- 
tionally, the  eclipse  obscuration  function  for  each  source  must  be  de- 
fined as  apptoptlato.  The  obscuration  functions  may  differ  for  each 
source,  and  none  ot  them  mav  be  well-represented  bv  the  visible  ob- 
scuration function. 

Sources  ot  ionization  f or  the  1'lftf.  eclipsed  D-region  are:  galactic 
cosmic  tavs,  solat  lvman-alpha  anti  Lvman-beta,  solar  x-ravs,  and 
energetic  electrons  precipitated  from  fix'  radiation  belts.  The  effects 
ot  photoloni zat ion  ot  0,('.\^)  bv  solat  ultraviolet  (UV)  radiation  were 

also  considered,  but  this  source  is  negligible.  The  Ionizing  sources 
for  the  l'*bt>  eclipse  utilized  in  this  report  were  obtained  from  measure- 
ments, for  example,  of  solat  x-tav  and  Lvman-alpha  fluxes,  and  from 
models.  The  sources  and  their  related  data  are  summarized  In  table  3. 

Cosmic  Rav  Ionization 

Galactic  cosmic  ravs  (CCK)  are  an  almost  constant  temporal  source  of 
ionization  in  the  upper  atmosphere.  Their  intensity  at  the  earth 
varies  with  geomagnetic  latitude  and  has  a small  temporal  modulation 
imposed  by  the  average  value  of  solar  wind  intensity,  which  is  con- 
trolled bv  solar  activity.  because  of  their  verv  high  energv,  their 
stopping  range  in  the  atmosphere  is  in  the  lower  stratosphere,  and 
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Figure  i',  Comparison  of  partial  reflection  data  with  in— situ  rocket  probe 
i'  lei- iron  density  measurements. 
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hence  their  ionization  profile  Is  proportional  to  atmospheric  neutral 
particle  density  at  higher  altitudes.  Heaps  [13]  outlines  the  deriva- 
tion of  the  Ionospheric  ionization  rate  due  to  GCR  and  presents  an 
empirical  formula  which  describes  the  geomagnetic  latitude  and  altitude 
variation  of  the  ionization  source  strength. 

The  GCR  ionization  source  intensities  below  60  degrees  magnetic  lati- 
tude at  solar  maximum  and  at  solar  minimum  are  described  by: 

4 -20 

Q(GCR:  S .MAX)  - (174  + 1930  sin  A)  x 10  N 
Q(GCR:  S .MIN)  - (174  + 2840  sin14  A)  x 10~:°  N 

_3 

where  A is  the  geomagnetic  latitude,  and  N is  the  number  density  (.cm  ) 
of  neutral  particles  in  the  atmosphere.  To  find  the  GCR  ionization 
source  for  the  circumstances  of  the  1966  solar  eclipse,  an  interpola- 
tion between  these  bounding  values  was  made  ns  described  in  Heaps  [13]. 
For  1966,  at  a geomagnetic  latitude  of  24  degrees, 

Q(GCR)  • RN 

— 1 8 — I 

where  K “ 2.32  x 10  ions  s . The  altitude  profile  of  the  cosmic 
ray  ionization  source  is  illustrated  in  fig.  7. 

Precipitating  Electrons 

A number  of  experimental  observations  and  theoretical  factors  suggest 
that  precipitating  electrons  introduced  a temporally  constant  source  of 
ionization  in  the  0-region  during  the  eclipse.  Evidence  of  trapped  or 
precipitating  electrons  was  obtained  by  rocketborne  x-ray  counters 
which  observed  anomalously  high  counting  rates  from  about  160  km  alti- 
tude to  near  apogee  both  before  and  during  the  eclipse  [14],  Absorp- 
tion measurements  made  during  totality  suggest  that  the  D-region  elec- 
tron densitv  is  anomalously  high  at  this  time,  which  could  only  result 
from  either  an  unknown  source  of  detachment  of  electrons  from  negative 
ions  or  from  an  added  ionization  source  [15],  Balloonborne  x-ray  de- 
tectors flown  in  this  region  have  detected  bremsstrahlung  from  pre- 
cipitated energetic  electrons  [16].  In  addition,  Potemra  and  Zmuda 
[l/]  showed  that  electron  prec ip i tat  ion  at  midlatitudes  could  contrib- 
ute a significant  nighttime  ionization  source  to  the  D-region,  In 
view  of  the  eclipse  experiment  location  on  the  western  portion  of  the 
South  Atlantic  geomagnetic  anomaly,  the  possibility  of  electron  pre- 
cipitation from  the  lower  trapped  radiation  belt  must  be  accounted  for. 

Eather  and  O'Brien  [18]  attempted  to  detect  electron  precipitation  in 
this  geographical  region  in  1967  by  photometric  methods  but  did  not 
detect  any  4278-angstrom  band  emission  enhancements  caused  by  electron 

ionization  and  excitation,  N-,  + e N^*  + 2e,  within  their  instrumen- 
tation sensitivity  limits  of  0.2  Rayleighs.  From  these  observations, 
they  concluded  that  the  particle  precipitation  over  the  region 
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Figure  7.  Altitude  dependent e of  the  ionizing  sources  for  the  12  November  1966  solar  eclipse 
pre-eclipse  conditions 
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accessible  to  their  view  was  less  than  7 x 10  4 ergs  cm  “s~  , 

4 -•>  -1 

10  electrons  cm  's  of  50  keV  energy.  As  Potemra  and  Zmuda  point  out, 
electron  fluxes  much  lower  than  tills  energy  limit  are  capable  of 
creating  D-reglon  ionization  which  is  measurable  in  the  absence  of 
normal  daytime  ionization  processes. 

In  the  absence  of  any  evidence  of  a temporally  unique  particle  precipi- 
tation event  occurring  during  the  eclipse  period,  the  precipitating 
electron  flux  was  computed  from  the  measured  trapped  electron  population 
on  the  L-shell  which  intersects  the  100  km  altitude  level  in  the  geo- 
graphical location  of  the  eclipse  measurements.  The  L-value  at  this 
location  is  about  1.2  which  corresponds  to  an  invariant  geomagnetic 
latitude  of  24  degrees.  The  average  electron  lifetime  for  this  L-shell, 
which  is  independent  of  energy,  is  250  days  [19]. 

The  trapped  electron  flux  in  several  energy  ranges  from  50  keV  to  1.7 
MeV  was  measured  by  Pfitzer  and  Winckler  [20]  on  30  September  1966. 

Given  the  average  loss  rate  from  this  time  to  the  eclipse  date,  with  no 
additional  replenishment,  about  88  percent  of  the  observed  flux  would 
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remain  and  about  4.1  x 10”  of  this  flux  would  be  lost  each  second. 

This  fraction  represents  the  minimum  average  precipitated  flux  over  the 
entire  longitude  range  of  the  L-shell  and  hence  represents  a minimum 
value  for  the  South  Atlantic  anomaly  region, 

Gassman  and  Pike  [21]  showed  that  the  eclipse  location,  on  the  westward 
side  of  the  anomaly,  can  receive  precipitating  electrons  from  over  300 
longitudinal  degrees  of  the  trapping  volume.  Because  the  mirror  point 
altitudes  for  a given  magnetic  moment  of  the  trapped  electrons  are 
lowest  and  are  descending  in  altitude  in  the  western  portion  of  the 
anomaly,  the  entire  average  electron  loss  rate  may  be  attributable  to 
loss  in  this  region  in  the  southern  hemisphere;  hence,  the  precipitating 
particle  fluxes  may  be  an  order  of  magnitude  greater  than  the  average 
value  would  indicate.  The  measured  trapped  fluxes  and  the  two  limiting 
cases  for  precipitating  fluxes  are  summarized  in  table  4. 

The  altitude  distribution  of  the  ionization  depends  to  a large  extent 
upon  the  energy  spectrum  of  the  precipitating  particles.  The  energy 
spectrum  as  measured  by  Pfitzer  and  Winckler  is  considerably  harder 
than  those  sample  spectra  adopted  by  Potemra  and  Zmuda  for  midlatitudes. 
For  example,  the  integral  energy  spectrum  derived  from  the  measurements 

is  approximately  F ■ E~*;  whereas  Totemra  and  Zmuda  adopt  energy 

dependences  of  L~ ' and  E-1.  However,  the  trapped  radiation  belts  do 
exhibit  harder  spectra  for  the  electron  flux  at  lower  L values  [22]; 
therefore,  the  measured  spectrum  is  used  to  compute  the  ionization 
prof ile. 

Ionization  profiles  for  both  cases,  uniform  loss  averaged  over  longi- 
tude and  loss  averaged  within  the  westernmost  60  degrees  of  the  South 
Atlantic  anomaly,  are  presented  in  fig.  7.  These  curves  were  computed 
following  the  method  of  Rees  [23]. 
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Solar  l.vman-A  1 pin  and  l.vman-Beta 


The  solar  l.vman-alpha  omission  intonsitv  is  a function  ot  solar  aotivitv 
( J | . Moasuromonts  ot  tho  solar  l.yman-a  l pha  energy  flux  .it  t imos  noar 

tin1  oolipso  range  t rom  1.  IS  org/om"  soo  to  S.7  ergs/em“  soo  |JS,Jti|. 

I'ho  nppor  valno  is  tnoro  oonsistont  with  that  expected  tor  tin'  solar 

aotivitv  lovol  at  the  t inn  ot  tho  oolipso.  lhe  lower  value  mav  have 

boon  produced  by  ionization  chamber  window  degradation  before  or  during 

> 

tho  measurement  (JS).  For  purposes  ot  this  work,  tho  a. 7 erg/cm"  soo 
value  is  adopted,  whieh  oorrosponds  to  a photon  t lux  ot  1.0  \ 10^' 
photons  otn"'  s”  ' . 

rho  solar  l.vman-hi'ta  flux  is  also  dependent  upon  tin1  level  ot  solar 
aotivitv.  Swider  |J7|  adopts  a value  ot  100  for  the  ratio  ot  l.vman- 
alpha  to  l.vman-beta  photon  fluxes,  whereas  Tousey  [JS]  reports  a rat  io 
ot  S>  t or  thi'  respeot ive  energy  l luxes.  No  direet  measurements  ot 
l.yman-beta  were  made  during  the  eo 1 ipse  period.  In  addition  to  direot 
l.yman-beta,  the  eontinuum  near  the  line  also  produoes  a small  degree  ot 
ionisation.  riierefore,  the  equivalent  l.vman-beta  tlux  is  taken  as 

1.0  x It)'1  photons  om”“  s-'. 

rho  ionization  produeod  by  solar  l.vman-alpha  and  l.vman-bt  ta  radiation 
depends  upon  the  altitude  distribution  of  the  Nit  and  0,  whieh  are  the 

ionized  speeies,  respeet  ive Iv . l.vman-alpha  is  strongly  absorbed  by  0,; 

lienee  the  distributions  of  both  neutral  speeies  at  feet  the  ionizing 
souree  strength  for  l.vman-alpha. 

Measurements  and  t heoret  leal  estimates  of  NO  density  versus  altitude 
have  varied  bv  more  than  a faetor  of  10.  Original  eomput.it  ions  ot  the 
ionosphorie  response  to  the  ll>Ph  solar  eel  ipse  1 1]  utilized  NO  profiles 
derived  t rom  Me i ra  (J9).  Measurements  bv  Baker  et  al.  | 10)  made  at 
White  Sands,  whieh  is  at  a corresponding  latitude  in  the  northern  hemi- 
sphere, indieate  an  overall  lower  NO  density  profile  with  a more  pro- 
nounced minimum  in  density  near  SS  km.  A comparison  of  the  ionizing 
souree  strength  utilizing  these  two  NO  profiles  is  presented  in  tig.  7, 
for  preeel  ipse  conditions.  The  figure  shows  that  for  the  first  souree 
strength  profile  (high  NO  content)  the  l.vman-alpha  ionization  predomi- 
nates over  a much  larger  altitude  range  than  for  the  low  NO  content 
eomput at  ion . 

Because  ot  the  wide  differences  in  empirical  and  theoretical  NO  pro- 
tiles, and  their  potential  ot foots  upon  0-region  ionization  profiles, 
the  single  "proper"  or  optimum  profile  as  an  a priori  input  to  our 
study  cannot  be  arbitrarily  selected.  Rather,  the  available  data  on 
electron  and  ion  density  profiles  have  been  used  and  an  attempt  has 
been  made  to  select  the  optimum  NO  profiles  based  upon  an  empirical 
best  fit  to  tlie  ionization  data.  This  procedure  is  outlined  in  the 
section  on  comparison  of  the  DA1RCHKM  code  with  the  19bb  eclipse  data  base. 
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The  solat  l.vman-beta  ionization  profile  is  computed  in  a straightforward 
manner  in  the  I'AIKi'Ui’M  code  and  depends  enlv  upon  the  0,  neutral  density 

prdtlle  and  the  overall  l \ can-bet  a flux.  The  ionization  source  strength 
versus  altitude  t ror.i  solar  Lvcian-beta  is  Illustrated  in  fig.  ?. 

Solar  X-Kays 

Measurements  o!  the  unobscured  solar  x-ra\  t lux  weie  made  Iron  satel- 
lites and  from  rockets  launched  i i ora  the  Brazilian  eclipse  site.  These 
measurements  nav  be  divided  into  two  spectral  groups  according  to  the 
legion  of  the  ionosphere  which  thev  ionize:  the  hard  x-iavs  between 
1 and  h'  angstroms  aie  capable  ot  ionizing  the  0-region  below  ll0  km; 
the  soft  x-ravs  between  about  10  and  100  angstroms  contribute  ioniza- 
t ion  mainlv  to  the  l -region  above  °0  km.  both  groups  are  important 
around  °0  to  100  km. 

The  measured  unobscured  x-ra\  t luxes  in  both  groups  are  summarised  in 
table  >.  N.'te  that  the  experimental  techniques  for  determining  these 
t luxes  depend  upon  the  x-rav  absorption  edges  ot  the  photometer  windows 
which  general lv  do  not  coincide  with  the  t wo  bands  specified.  There- 
fote,  the  spectral  distribution  and  total  intensity  within  these  bands 
must  be  inf ei red  Item  the  data. 

The  spectral  distributions  in  the  two  x-rav  specttal  regions  derived 
bv  Swidei  (O'!  wete  adopted  toi  tin1  present  studv.  These  spectral 
distributions  were  normalized  to  correspond  with  the  mote  limited  band- 
pass x-ru\  t luxes  measured.  In  general,  details  ot  the  spectra  derived 
bv  Swider  and  the  measurements  for  the  eclipse  period  were  in  good 
agreement,  and  normal izat ion  ot  the  two  derived  spectra  to  the  measured 
x-rav  fluxes  should  be  accurate  to  within  about  JO  percent.  Figure  8 
compares  the  observed  fluxes  and  the  x-iav  spectra  in  the  two  wave- 
length regions  adopted  tot  this  studv. 

The  Ionization  profile  resulting  t rom  absorpt ion  ot  the  t wo  solar  x-rav 
emission  bands,  l to  10  angstroms  and  SI  to  100  angstroms,  was  computed 
tot  the  unobseured  situ.  After  normal i za t ion  foi  the  incident  integrated 
unobscured  solar  t lux  in  the  two  bands  of  interest,  the  ionization  pro- 
files computed  were  net  significantly  ditterent  from  those  computed  bv 
Swidei  toi  a solar  zenith  angle  of  JS  degrees.  The  computed  x-rav 
ionization  source  strength  profiles  are  Illustrated  in  tig. 

Pet  tuition  ot  the  Kcllpse  Obscuration  Functions 

The  temporal  variation,  in  addition  to  the  absolute  magnitude  versus 
altitude,  of  the  ionizing  radiations  must  be  defined  during  the 
eclipse.  Solar  radiation  eclipse  obscuration  functions  were  defined 
for  solar  Ivman-alpha  and  Lvman-beta,  visible  and  near  VY  radiation, 
and  both  the  hard  and  sot t solar  x-rav  emissions. 

Visible  and  tieai  TV  obscurations  functions  are  computed  in  PAlKl'HKM  bv 
using  the  geomet i teal  tormula  of  Rvdheck  [ M ' and  the  I’V  limb  darkening 
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These  obscuration  functions  .are 


model  of  Tierce  and  Waddell  [32]. 

I 1 lust  rated  in  I i|',.  9. 

L.Vman-alpha  and  l.vman-beta  obscuration  functions  are  somewhat  more  com- 
plex. before  1958,  the  i lux  ol  l.yman-alpha  radiation  Incident  upon  the 
atmosphere  within  the  region  of  totality  was  believed  to  be  negligible. 
On i Inc,  the  1 9M>  solar  eclipse,  two  independent  measurements  ol  l.yman- 
alpha  l lux  within  the  umbra  were  made  and  both  showed  a residual  flux. 
Measurements  made  hv  HI ament  and  Malique  ( 2 S ] showed  solar  limb 
bt  ightening  in  1. 1 which  extended  above  the  photosphere  or  optical  1 litih 

to  altitudes  greater  than  2500  km.  A residual  L(t  flux  ot  0.14  percent 

was  reported  lor  totality.  Bowhiil*  estimated  a residual  flux  of  about 
0.5  percent,  based  upon  dltterent  data. 

Two  possible  mechanisms  mav  account  for  residual  I,  flux  observed  at 

a 

totality:  resonant  scattering  lit  the  uneclipsed  and  penumbral  portions 

ot  the  lerrosti ial  hydrogen  geocorona,  and  failure  to  achieve  complete 
totality  over  the  I.u  disk.  Sears  | 1]  estimated  the  magnitudes  of  these 

respective  components  as  being  loss  than  0.5  percent  and  about  0.1  per- 
cent, respectively.  for  litis  study,  the  residual  La  I lux  at  totality 

as  0.5  percent  was  used.  The  temporal  behavior  of  the  L ionizing 
source  strength  during  the  eclipse  Is  summarized  in  lig.  9. 

Treatment  of  the  l.ynian-heta  obscuration  function  is  considerably  simpler 
titan  for  I.vman-alpha  because  the  hydrogen  geocorona  is  much  less  opti- 
cally thick.  Motor  ( 13]  has  shown  that  for  near  zero  solar  zenith 
angle  and  lor  exospheric  temperatures  above  about  900  °K,  the  optical 
thickness  through  the  geocorona  to  an  altitude  of  100  km  is  less  than 
two.  Honce,  the  requirements  for  a radiation  transport  model  are  con- 
siderably lessened.  The  radiation  transport  considerations  for  this 
obscuration  fund  Lon  are  neglected  and  a visible  obscuration  function 
as  indicated  in  t ig.  4 will  provide  results  of  adequate  accuracy. 

X-Kav  Obscuration  Function 

The  x-ray  obscuration  function  depends  stronglv  upon  the  localized  dis- 
tribution of  x-ray  sources  on  the  solar  disc.  An  analytical  function 
cannot  be  derived  in  general,  and  experimental  data  must  be  sought. 
Photographs  of  the  solar  disc  in  x-ray  wavelengths  were  obtained  by 
Ihulerwood  et  al . |34]  for  the  unohscured  sun  at  the  time  ol  the  solar 
eclipse.  Although  it  is  well-known  and  strongly  evident  in  the  x-ray 
photographs  that  the  harder  portions  of  the  x-ray  spectrum  arc  more 
localized  than  the  soft  emission  spectrum,  thus  producing  a larger 
temporal  gradient  i n ionization  source  strength,  this  analysis  shows 
that  such  temporal  gradients  cannot  be  resolved  adequately  by  the  U- 
regioti  measurement  techniques  and  the  available  data.  Hence,  although 


^Comments  made  at  the  1966  Solar  Kclipse  Symposium,  Sao  .lose  dos 
Campos,  Brazil,  February  1968. 
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not  justified  in  detail,  the  soft  x-ray  obscuration  function  derived 
from  the  27-  to  40-angstrom  photograph  is  used  for  the  entire  soft 
x-ray  flux. 

Measurements  of  solar  x-rays  during  the  1970  eclipse  by  Accardo  et  al. 
[35]  give  residual  values  of  16  percent  during  totality  for  the  x-ray 
band  from  44  to  60  angstroms,  which  is  in  very  good  agreement  with  the 
obscuration  function  used  here.  Measurements  in  the  2-  to  8-angstrom 
band  showed  a residual  value  of  5 percent  during  totality.  Subse- 
quently, the  hard  x-ray  obscuration  function  used  here  follows  the 
soft  x-ray  function  down  to  about  5 minutes  before  totality  and  then 
descends  to  the  lower  5 percent  value  during  totality. 

Summary  of  Ionizing  Source  Strength  and 
Obscuration  Function  Calculations 

One  may  conclude  from  the  results  of  the  definition  of  the  ionizing 
sources,  their  altitude  variations,  and  eclipse  obscuration  functions 
that  an  adequate  description  of  the  ionizing  source  input  to  the  D- 
region  under  eclipse  conditions  entails  a complicated  superposition  of 
time  and  altitude-varying  functions.  It  is  instructive  to  compare  the 
altitude  profiles  of  ionizing  source  strength  during  totality  in  a for- 
mat equivalent  to  that  useu  in  fig.  7 for  preeclipse  conditions. 

Figure  10  illustrates  these  results.  The  figure  illustrates  the  im- 
portance of  the  precipitating  electron  flux  over  the  entire  lower  D- 
region  compared  with  the  eclipsed  solar  ionizing  sources  such  as 

and  x-rays.  As  will  be  shown  later,  a significant  precipitating  elec- 
tron flux  is  required  to  model  the  ionization  profile  in  preeclipse 
conditions.  As  in  the  preeclipse  condition,  the  galactic  cosmic  ray 
ionization  is  negligible  over  the  altitude  range  of  interest  to  this 
study. 

DESCRIPTION  OF  THE  DA1RCHEM  COMPUTER  CODE 


DAIRCHEM  determines  the  number  density  as  a function  of  time  at  a 
specified  altitude  for  each  of  several  charged  and  neutral  species.  A 
large  set  of  time-dependent,  coupled  chemical  reaction  equations  is 
solved,  of  the  form 


q^  + 1'  (rates  forming  y^)  - 1'  (rates  removing  y^) 


where  y.  is  the  time-dependent  number  density  of  the  i*"*1  species; 

is  the  nonchemical,  nonphoton  source  term;  and  the  last  two  terms 
represent  formation  and  removal  terms  due  to  chemical  and  photon  re- 
actions. The  current  reaction  set  contains  493  reactions,  involving 
64  neutral,  positive,  and  negative  species. 

Basic  routines  for  solving  these  stiff,  coupled  differential  equations 
have  been  taken  from  the  existing  AIRCHEM  atmospheric  chemistry  code 
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Figure  10.  Altitude  dependence  of  ionizing  sources  for  the 


[36].  The  latest  version  of  the  numerical  integration  routine  is  based 
on  a predictor-corrector  algorithm  called  the  K-method  137,38],  Its 
Improved  speed  can  now  simulate  a 24-hour  period  at  any  mesospheric 
altitude  In  less  than  2 minutes  CPU  time  on  a CDC  7 600  computer. 

The  solar  1 lux  arriving  at  the  altitude  of  interest  Is  determined  by 
taking  into  account  absorption  by  N., , COj,  and  Oj.  The  column  den- 
sity ot  each  absorbing  species  is  calculated  along  the  geometric 
straight-line  path  traveled.  Solar  zenith  angles  near  and  greater  than 
4(1  degrees  are  readily  accommodated.  Scattering  into  the  volume  ol 
interest  and  away  from  the  line  of  sight  are  not  presently  incorporated. 
The  absorption  ot  solar  radiation  and  dissociation  of  Ot  in  the  175  to 

210  nm  spectral  region  dominated  by  the  Schumann-Runge  bands  are 
handled  by  using  the  semlempirical  method  of  Kockarts  [39].  When  a 
solar  eclipse  is  in  progress,  the  incident  solar  flux  is  (list  reduced 
by  the  appropriate  obscuration  functions  for  the  visible,  UV,  l.vman- 
alpha,  and  soft  and  hard  x-ray  port  ions  of  the  spectrum. 

Photodetachment  ol  electrons  from  negative  ions,  and  also  photo- 
dissociation  of  ions  and  neutrals  are  considered,  and  the  appropriate 
cross  sections  have  been  taken  from  the  literature  and  sponsored  re- 
search work  (e.g. , Vamlerhot  I [40],  Vanderhot'l  and  Rover  [411,  Cosby  et 
al.  [42]  and  references  contained  therein).  In  practice,  chemical 
means  of  detaching  electrons  from  negative  ions  are  found  to  dominate 
the  photodetachment  processes  [43]. 

A more  detailed  description  ot  the  PA1RCHKM  computer  code  as  it  is  used 
in  the  current  set  of  runs  mav  be  found  in  Uoock  and  Heaps  [44]. 

COMPARISONS  OF  DAIRCHEM  WITH  THE  196b  ECLIPSE 
DATA  RASE  PRKKCL1PSE  CONDITIONS 

DAIRCHEM  runs  were  made  at  altitudes  between  60  and  9(1  km  to  develop 
profiles  of  electron  and  total  positive  ion  densities  to  compare  with 
experimental  measurements.  To  accurately  simulate  the  altitude  pro- 
files of  electron  densities  during  uneclipsed  conditions,  the  precipi- 
tating electron  ionization  source  previously  described  had  to  be  in- 
cluded. The  magnitude  ot  this  source  of  ionization  was  determined  bv 
comparing  code  results  at  60  and  66  km  with  measured  electron  densities. 
The  results  of  this  procedure  showed  that  the  magnitude  of  the  precipi- 
tating electron  ionization  source  was  the  maximum  value  modeled  (fig.  7). 
Above  66  km,  the  solar  l.yman-alpha  ionizing  source  becomes  increasingly 
important,  depending  in  magnitude  and  profile  upon  the  NO  model  chosen. 
The  effects  of  varying  the  Nil  density  profile  modes  were  illustrated  in 
fig.  7.  Rest  fit  ot  the  code  predictions  to  the  measured  electron 
densities  over  the  70  to  80  km  altitude  range  occurred  when  the  Raker 
et  al.  [29]  profile  was  chosen,  although  a fairly  wide  variation  of  NO 
densities  in  this  region  can  be  accommodated . Figure  11  illustrates 
the  simulated  and  experimentally  determined  electron  and  total  positive 
ion  altitude  profiles  for  preeclipse  conditions. 


ELECTRON  AND  TOTAL  POSITIVE  ION  NUMBER  DENSITY  (CM"3) 

Figure  11.  Electron  density  and  total  positive  ion  density  during  pre-eclipse  conditions  are  co-.pared 
for  DAIRCHEJi  predictions  and  rocket  probe  data. 
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While  the  electron  density  profiles  are  in  reasonable  agreement  over  the 
entire  altitude  range,  the  code  predicts  substantially  lower  values  of 
total  positive  ion  density  than  were  observed  over  the  60  and  83  km 
interval.  The  figure  shows  that  significant  negative  ion  densities  may 
be  inferred  as  high  as  90  km  from  the  data,  whereas  the  code  predicts 
very  low  negative  ion  densities  above  about  70  km.  A somewhat  different 
picture  emerges  if  the  total  positive  ion  density  curves  of  Narcisi  et 
al.  [10]  are  used  (figs.  3a  and  3b).  The  figures  show  that  the 
uneclipsed  positive  ion  density  remains  almost  constant  at  about  1 x 10^ 

cm  between  72  and  84  km  and  follows  the  electron  density  curve  above 
82  km.  These  data  imply  a negative  ion  cutoff  at  82  km;  however,  the 
accuracy  of  the  total  positive  ion  profile  below  about  90  km,  which  is 
inferred  from  the  ion  mass  spectrometer  measurements,  is  subject  to 
some  uncertainty.  Hence,  a quantitative  discrepancy  between  code  simu- 
lation and  observational  data  remains  with  respect  to  total  positive  ion 
density  and  inferred  negative  ion  density  in  the  60  to  85  km  region. 

Identification  of  Positive  Ion  Species 

The  positive  ion  species  measured  by  mass  spectrometer  probes  [10]  be- 
fore and  during  the  eclipse  may  be  compared  with  DAIRCHEM  predictions. 
Figures  5a  and  5b  illustrate  the  measured  positive  ion  species  profiles 
for  a certification  round,  several  days  prior  to  the  eclipse  and  at 
totality.  The  figures  show  that,  without  regard  to  the  magnitude  of 

the  total  positive  ion  density,  at  higher  altitudes  the  N0+  and  0.,+ 

species  predominate  throughout  the  eclipse,  whereas  below  about  82  to 
84  km,  hydrated  species  predominate. 

DAIRCHEM  simulations  of  positive  ion  densities  for  preeclipse  and  total- 
ity conditions,  shown  in  figs.  12a  and  12b,  correctly  give  the  observed 

N0+  and  0.,+  predominance  at  higher  altitudes.  The  simulated  0-,+  compo- 
nent is  somewhat  suppressed,  indicating  that  the  hard  x-ray  flux  or 
flux  of  precipitating  electrons  of  approximately  40  keV  energy  may  have 

been  underestimated.  DAIRCHEM  predicts  that  the  first  hydrate  of  N0+ 
is  a major  ion  in  the  80  to  86  km  region,  particularly  during  totality, 

but  N0+*H-,0  was  detected  on  only  one  flight  made  at  80  percent  obscura- 
tion (not  shown).  Two  possibilities  for  the  discrepancy  are  that  the 
rate  along  the  hydration  chain  is  too  high  (temperature  and  water  vapor 
content  were  not  measured),  and  that  the  hydrated  ions  were  broken  up 
by  passage  through  the  probe  shock  wave  or  by  wall  effects,  thus 

appearing  as  N0+.  Better  qualitative  agreement  is  obtained  between 
measured  and  simulated  N0+  profiles  if  the  simulated  N0+  and  N0+*Hn0 
are  added  together. 

At  altitudes  below  the  N0+  region,  hydronium  hydrates,  H^0+‘(H10)^, 

predominate  both  before  and  during  the  eclipse.  The  data  show  that 
masses  37  and  55  (n  • 1 and  2)  predominate,  whereas  DAIRCHEM  predicts 
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that  higher  hydrates  are  relatively  more  important.  Moreover,  o;  total- 
ity is  approached,  the  lower  hydrates  tend  to  grow  into  the  ligh«r  hy- 
drates; the  total  ion  concentration  remains  reasonably  constant.  T 1 * 
authors  believt  that  the  DAIRl'HKM  predictions  are  consistent  will,  the 
data  insofar  as  the  rapid  passage  of  the  ion  mass  spectrometer  payload 
through  the  lower  D-region,  at  supersonic  velocities,  will  tend  i.  dis- 
sociate the  weakly  bound  multiple  hydrated  species.  Hence,  h\dr..tid 
ions  of  mass  55,  73,  or  greater  tn.iv  be  observed  as  lower  hydrates,  -5, 
37,  or  1H. 

The  major  d i t f >.  renves  then  between  the  simulated  and  measured  positive 
ion  densities  are  the  quantitative  failure  ot  the  code  to  yield  the 
large  total  positive  ion  densities  which  were  inferred  from  rocket 
probe  moasut ement s and  the  tendency  of  the  code  to  predict  highet  hy- 
drates of  N0+  and  11  than  were  measured. 

Temporal  Variations  ot  Flection  Density 

One  of  the  major  reasons  for  undertaking  the  detailed  comparison  ot 
DAIRCHEM  with  the  l^tib  solar  eclipse  experiment  is  to  model  the  very 
rapid  changes  of  electron  density  observed  to  occui  immediate! y at 
second  contact.  Previous  ground-based  measurements  reported  by  Sears 
[7,45]  had  indicated  that  there  were  rapid  and  large  changes  in  D- 
region  electron  density  associated  with  the  totality  region  which 
appeared  to  be  greater  in  magnitude  than  could  be  easily  explained  by 
lumped  parameter  reaction  rate  codes. 

The  temporal  changes  In  electron  density  measured  by  the  several  rocket 
flights  within  a few  hundred  seconds  of  totality  were  illustrated  in 
fig.  3.  The  figure  shows  a reduction  of  electron  densitv  et  a factor 
of  3 to  10  within  a few  tens  of  seconds  or  less  at  second  contact.  The 
figure  also  shows  a clear-cut  and  rapid  decline  in  electron  density  at 
second  contact  based  not  only  upon  the  combined  data  but  also  upon  com- 
parison of  measurements  made  within  individual  groups.  Further,  the 
cross-calibration  of  rocket  probe  results  with  each  other,  and  with 
ground-based  electron  density  measurements,  shows  that  the  observed 
density  decrease  is  much  greater  than  the  experimental  errors  inherent 
in  the  Individual  measurement  technique. 

The  attempt  of  the  DAIRCHEM  code  to  adequately  simulate  the  rapidity 
and  magnitude  of  the  electron  density  decrease  at  second  contact  Is 
illustrated  in  fig.  13.  The  Inability  of  the  code  to  produce  any  rapid 
drop  In  electron  density  above  70  km  during  totality  is  apparent. 

To  determine  the  reason  for  the  computer  code's  failure  to  simulate  the 
rapid  changes  In  electron  density  about  totality,  a closer  look  at  the 
negative  ion  chemistry  must  be  taken.  A schematic  chart  of  the  flow  of 
charge  between  electrons  and  negative  Ions  la  Illustrated  In  fig.  14, 
with  a more  detailed  listing  of  pertinent  reactions  given  In  table  6. 

The  moat  rapid  electron  loss  process  is  the  three-body  attachment  to  CU. 
After  Oj-  Is  formed,  further  reactions  carry  the  negative  charge  down 
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TABLE  6.  IMPORTANT  ELECTRON  ATTACHMENT /DETACHMENT 
AND  CHARGE  TRANSFER  REACTIONS 


Reaction 

Critical 

Minor  Species 

Rate  Constant 

0 j + e -*•  O”  + Oj 

°3 

9.0 

(-12)(T/300)1'5* 

0.,  + 0,,  + e -►  0-,  + 0^ 

- 

1.4 

(-29)(T/300)_1  exp(-600/T) 

02  + + e -►  0^“  + N2 

- 

1.0 

(-31) 

o2  + o -♦  o”  + o2 

0 

1.5 

(-10) 

o’  + 0 + 0 ' + 0 

3 3 

°3 

5.3 

(-10) 

0,"  + 03  -►  o3"  + 0., 

°3 

6.3 

(-10) 

03"  + co.,  -►  C03"  + 0, 

- 

5.5 

(-10) 

co3~  + 0 -•>  0o~  + co,, 

0 

1.1 

(-10) 

0 + 0 -♦  0^  + e 

0 

2.0 

(-10) 

0+0  (*A  ) -*■  0,  + e 

2 8 3 

0,  (*A  ) 

2 8 

3.0 

(-10) 

0^  + 0 -»  03  + e 

0 

1.5 

(-10) 

°2~  + 02  (*Ag)  ->  20,  ♦ e 

0,  (*A  ) 

2 8 

2.0 

(-10) 

♦Read  as  9.0  x 10  (T/300)1'^;  general  form  A (T/300)8  exr  (C/T) 
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the  chain  (principally  through  the  intermediate  ion  0^  to  CO  ) to  the 
more  stable  negative  ion  NO^-,  its  hydrates  and  clusters.  No  rapid, 

direct  attachment  processes  of  electrons  to  the  more  stable  negative 
ions  are  currently  known. 

The  rate  of  electron  three-bodv  attachment  is  rapid  enough  to  explain 
the  type  of  electron  density  fluctuations  shown  in  fig.  13,  but  the 
rate  of  electron  detachment  is  equally  as  fast.  The  species  responsible 

for  detaching  electrons  from  0^'  and  0~  are  0 and  0o(1A  ).  If  there 

*■  R 

were  rapid  decreases  in  these  latter  densities  at  totality,  then  per- 
haps the  electrons  could  stay  attached  to  the  0,”  and  0-,  thus  account- 
ing for  the  rapid  decrease  in  the  free  electron  density.  Similarly, 
the  crucial  minor  neutral  species  responsible  for  transfer  of  charge 

from  0,"  to  C03_  is  0^.  If  Oj  were  to  rapidly  increase  during  totality, 

then  charge  could  rapidly  pass  down  the  negative  ion  chain  and  no  longer 
be  affected  by  the  rapid  detachment  processes.  However,  the  time  con- 
stants for  the  minor  neutral  constituents  are  generally  longer  than 

those  for  ionic  species,  so  that  the  changes  in  the  0,  and  0^ 

densities,  shown  in  fig.  15,  are  relatively  smooth  and  only  ajnount  to 
factors  of  2 or  3.  The  variations  needed  to  explain  the  observed  elec- 
tron density  fluctuations  would  have  to  be  order  of  magnitude  changes 
on  a time  scale  of  100  seconds. 

Variations  in  the  known  reaction  rates  to  simulate  a rapid  equivalent 
two-body  attachment  process  involving  known  species  have  also  been  tried 
[46],  but  with  unsatisfactory  results  in  simulating  both  the  rapidity 
and  magnitude  of  the  electron  density  decreases. 

Therefore,  no  appeal  can  be  made  to  the  known  gas-phase  chemistry, 
which  controls  the  electron  attachment  and  detachment  processes,  to 
adequately  explain  the  rapid  changes  in  electron  density  in  the  65  to 
85  km  region. 

SUMMARY  AND  CONCLUSION 

The  results  of  comparing  the  predictive  capabilities  of  the  DAIRCREM 
code  with  the  data  obtained  on  the  1966  solar  eclipse  are  summarized  as 
follows : 

1.  The  code  provides  a good  fit  to  the  electron  density  versus 
altitude  for  preeclipse  undisturbed  conditions. 

2.  The  code  prediction  of  the  total  positive  ion  density  in  the 
altitude  range  from  about  60  to  85  km  at  all  times  is  too  low  by  a 
factor  of  2 to  5. 

3.  The  dynamical  reduction  of  the  electron  density  at  second  con- 
tact at  altitudes  below  about  85  km  is  not  predicted  by  the  code. 
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Figure  15.  Variation  during  the  eclipse  of  several  of  the  minor  neutral 
species  important  in  formation  of  negative  ions  (Oj)  and 
electron  detachment  (0  and  )). 


4.  There  are  physical  limitations  on  adjusting  the  reaction  rates 
and  processes  within  the  code  such  that  the  eclipse  data  can  be  fit, 
especially  near  second  contact.  Since  the  current  gas  phase  chemistry 
cannot  model  the  actual  events,  the  potential  surface  chemistry  effects 
of  ambient  mesospheric  aerosols  may  have  to  be  included. 

At  this  point  it  is  not  clear  whether  the  discrepancy  between  prediction 
and  measurement  with  respect  to  positive  ion  densities  or  of  temporal 
electron  density  behavior  in  the  t>()  to  90  km  altitude  range  is  a failure 
in  the  reaction  chemistry  contained  in  the  code,  or  is  caused  by  other 
physical  processes  in  the  ionosphere  which  are  not  modeled  at  all.  The 
very  rapid  fall  in  electron  density  at  second  contact  could  be  explained 
by  attachment  to  0-i,  but  only  if  the  subsequent  detachment  processes  are 

ignored.  Another  possibility  for  such  rapid  changes  is  by  loss  to 
another  species,  for  example,  by  surface  attachment  to  an  aerosol  spe- 
cies. The  large  difference  between  the  observed  electron  and  positive 
ion  density  in  the  b 0 to  90  km  region,  coupled  with  the  fact  that  the 
current  gas  phase  chemistry  predicts  verv  few  negative  ions  above  70  km, 
would  indicate  that  negatively  charged  aerosols  are  a distinct  possibil- 
ity. 

The  inclusion  of  aerosols  would  require  the  introduction  of  a major  new 
set  of  physical  and  chemical  processes  into  the  conventional  mesospheric 
physics  and  chemistry,  and  verification  by  both  field  and  laboratory 
experiments . 
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